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The adiabatic perturbation is shown to be singular from the exact solution of a spin-1/2 particle 
in a uniformly rotating magnetic field. Due to a non-adiabatic effect, its quantum trajectory on 
a Bloch sphere is a cycloid traced by a circle rolling along an adiabatic path. As the magnetic 
field rotates more and more slowly, the time-energy uncertainty, proportional to the distance of 
the quantum trajectory, calculated by the exact solution is entirely different from the one obtained 
by the adiabatic path traced by the instantaneous state. However, the non-adiabatic Aharonov- 
Anandan geometric phase, measured by the area enclosed by the exact path, approaches smoothly 
the adiabatic Berry phase, proportional to the area enclosed by the adiabatic path. The singular 
limit of the time-energy uncertainty and the regular limit of the geometric phase are associated 
with the arc length and arc area of the cycloid on a Bloch sphere, respectively. Prolate and curtate 
cycloids are also traced by different initial states outside and inside of the rolling circle, respectively. 

The axis trajectory of the rolling circle, parallel to the adiabatic path, is shown to be an example of 
transitionless driving. The non-adiabatic resonance is visualized by the number of complete cycloid 
arcs. 


Perturbation theorjsii^ is widely used in many fields 
of science and engineering as an effective method to find 
an approximate solution to a given problem, expressed 
in terms of a power series in a small parameter. In reg¬ 
ular perturbation calculations, one only keeps the first 
few terms of the expansion to obtain a good approxi¬ 
mate solution to the exact one, as the small parameter 
goes to zero. However, there are many interesting prob¬ 
lems that have no such uniform asymptotic expansion. 
These involve singular perturbations^^— . A classic exam¬ 
ple of this singular perturbation is the flow in the limit of 
zero viscosity^. When the viscosity, a small parameter, 
approaches zero, the solution of the Navier-Stokes equa¬ 
tion gives a completely different solution from the one 
obtained by taking the zero viscosity to start out with. 

Adiabatic perturbation^ is one of the fundamental ap¬ 
proximations used in many fields. Its classic applications 
include the Born-Oppenheimer approximation^ of decou¬ 
pling for the fast electronic motion from the slow ionic 
one, and adiabatic quantum computationiS, an alterna¬ 
tive to the quantum circuit model for quantum comput¬ 
ing. The adiabatic theorem dictates that as long as a 
system changes slowly enough, a quantum system start¬ 
ing from an eigenstate would remain in the instantaneous 
eigenstate of the time-dependent Hamiltonian up to the 
dynamical and Berry phase o^^i^^ . It may seem quite rea¬ 
sonable then that all physical properties in the adiabatic 
limit should be obtainable from the instantaneous eigen¬ 
state. However, this conjecture has never been proved. 

In this paper, we reveal singular features of the adia¬ 
batic approximation by studying the quantum dynamics 


of a spin-1/2 particle in a uniformly rotating magnetic 
field. Its quantum trajectory is shown to be a cycloid on 
the Bloch sphere, traced by a point on a rolling circle, 
of a radius determined by the angular speed of the mag¬ 
netic field, along the adiabatic path of the instantaneous 
eigenstate. We find the two basic geometric quantities, 
the distance and the enclosed area of the quantum trajec¬ 
tory, approach the different limits in the adiabatic limit. 
As the rotation of the magnetic field is slowed down, 
the non-adiabatic Aharonov-Anandan (AA) phasei^, the 
area enclosed by the quantum trajectory, goes to the 
adiabatic Berry phase, the area enclosed by the adia¬ 
batic path. However, the time-energy uncertainty, the 
distance of the quantum trajectory, does not converge to 
the minimum time-energy uncertainty of the adiabatic 
path. This singular feature of the adiabatic approxima¬ 
tion is explained by the arc length and arc area of the 
cycloid. In addition, the cycloid curve neatly explains 
some interesting physical results. First, the axis trajec¬ 
tory of a cycloid is interpreted as a transitionless driving 
that makes the quantum evolution follow the adiabatic 
path. Second, the non-adiabatic resonance condition is 
visualized by the number of perfect arcs of the cycloid. 
Finally, the exact cycloid, curtate and prolate cycloids on 
a Bloch sphere are generated by different initial states. 
Our results could be tested with a single qubit, a neutron, 
or light polarization, and could have important implica¬ 
tions for the application of the adiabatic perturbation, 
for example, adiabatic quantum computing and adiabatic 
quantum dynamics. 

Results 
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Spin-1/2 particle in a rotating magnetic field. We consider 
one of the simplest quantum systems, a spin-1/2 parti¬ 
cle in a rotating magnetic field B(t) = Bn{t) where we 
assume its strength B is constant and its direction n(t) 
rotates with constant angular speed w. Quantum dy¬ 
namics is governed by the time-dependent Schrodinger 
equation 


, ( 1 ) 

where the time-dependent Hamiltonian is given by the 
Syman interaction 


H{t) = 



( 2 ) 


Here cjq is the Larmor frequency (for electron spin tag = 
egB/2m), cr is the Pauli spin vector. As Feynman et 
al^ showed, with a Bloch vector r{t) = {tj}{t)\cr\tl;{t)) 
Equation o can be written as the dynamics of spinning 
top 


^ xr. (3) 

Before discussing the solution of Equation o, let us re¬ 
call adiabatic dynamics of a spin-1/2 particle. Adiabatic 
theorem dictates that when an applied magnetic field 
changes slowly enough, a quantum state \'ip{t)), adiabati- 
cally evolved from an initial instantaneous state |nj-(0)), 
remains in an instantaneous eigenstate |n±(t)) up to the 
dynamical phase and the adiabatic geometric phase 
7 ± (called the Berry phase) ^^d^ 

0 ::^ exp {±i^t + ij±) |n±(t)). (4) 

The instantaneous eigenstates |n±(t)) are the solution of 
H{t)\n±{t)) = =F^^|n±(<)) and written as |n+(t)) = 
cos||0) -I- e^^sinlll) and |n_(t)) = — sin||0) -b 

e“^ cos 111). Here 0 and (j) are the polar and azimuthal 
angles of n, respectively. In adiabatic limit, the Bloch 
vector r(t) = (i/)ad(t)|o’|^/’ad(t)), i.e., the spin direction is 
clearly align to the direction of the magnetic field n(t) 
in the parameter space. The Berry phase is expressed in 
terms of the geometric quantity, the solid angle S sub¬ 
tended by n(t) as 7 ± = 

The question we would like to explore is how the non- 
adiabatic trajectory of the Bloch vector r approaches to 
that of n when the magnetic field rotates slowly. To this 
end, the exact solution of Eq. o is obtained by trans¬ 
forming it to the adiabatic frame via the transformation 
|^/>(<)) = A{t)\(f{t)) where A{t) is composed of the col¬ 
umn vectors |n+(t)) and |n_(t)). In the adiabatic frame, 
the time-dependent Schrodinger Eq. becomes 

= Heslifiit)), (5) 


where the effective Hamiltonian is decomposed into the 
sum of the adiabatic and the non-adiabatic terms 

BA 

Hes = A-^HA - ihA-^ — 
dt 

hwo fl4> \ 

= -^0-2 -b —(I - cosOa^ + sinOa^) - —ay . 

( 6 ) 

For the magnetic field rotating with constant angu¬ 
lar speed, the effective Hamiltonian (jS]) becomes time- 
independent and Eq. o is exactly solvable. 

We consider the two cases: (i) the rotation along the 
latitude, (p = uit and constant 0, and (ii) the rotation 
along the meridian, 9 = tut and constant ip. Thus, in 
the adiabatic frame, the Bloch vector r{t) rotates with 
the frequency B around the new axis e at acute de¬ 
viated angle a relative to the z axis (adiabatic axis). 
The frequency B and the deviated angle a are given by 
B = y^Wq -b and a = tan“^ (w/wq) for the rotation 
along the meridian, and -b 2ujquj cos 6 + uj"^ and 

a = tan“^ [(tag -b w cos 0)/ta sin 0] for the rotation along 
the latitude, respectively. The slowness of the rotation of 
the magnetic field is relative to the Lamor frequency. So 
the ratio of two frequencies, A = uj/ujq controls the adi- 
abaticity of quantum dynamics. In limit of A —> 0, the 
quantum dynamics becomes adiabatic, and the radius a 
of the (imaginary) rolling circle along the adiabatic path, 
given by a = sin a, becomes smaller. The exact solution 
to Eq. o is written as 

IV'(t)) = A{t) A(0)|V'(0)). (7) 



Figure 1. Cycloids on a Bloch sphere. As the rotating mag¬ 
netic field traces the blue line, the Bloch vector makes various 
trajectories in red: (a) the exact cycloid, (b) prolate cycloid, 
(c) curtate cycloid, and (d) the axis trajectory, depending on 
the initial conditions. 

Cycloids on a Bloch sphere. We calculate the trajectory of 
the Bloch vector r{t) on the Bloch sphere with the exact 
solution Eq. ©• Fig. m plots the cycloids traced by the 
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Bloch vector r{t) = {ilj{t)\a\'tjj{t)) on a Bloch sphere for 
different initial states when the magnetic field is rotated 
clockwise with angular speed uj = 0.5 tag around the z- 
axis by azimuthal angle 9 = tt/S. On a plane, a cycloid is 
the curve traced by a point on the rim of a circle that rolls 
along a straight line. In classical mechanics, it is the solu¬ 
tion to two famous problems: brachistochrone (shortest- 
time) curves and tautochrone (equal-time) curvesi^*^. 
Like on plane, a cycloid on the Bloch sphere is traced 
by a point of an imaginary circle which rolls along the 
base line. Here the adiabatic path n(t), i.e., the trajec¬ 
tory of the magnetic field, plays a role of the base line, as 
shown by the blue curve in Fig. [T] The imaginary circle 
rolling along the adiabatic path represents non-adiabatic 
quantum dynamics. The radius of the rolling circle is 
determined by two frequencies wq and w, and given by 
a = sin a. Slower the rotation of the magnetic field, 
smaller the rolling circle. This illustrates clearly how the 
quantum trajectory approaches the adiabatic path in the 
adiabatic limit. 



plane cycloid 

spherical cycloid 

base line 

straight line 

adiabatic path n{t) 

circle radius 

a 

a = sin a 

rolling speed 

ifi = iit 

17 = \J uJq A- ixP' 

equation | 

' X = aip — b sin ip 
y — a — b sin ip 

7^=-^n.a|V>> 

f = -cao n X r 

curtate/prolate 

in/outside 

in/outside 

arc length 

8o 

4acosa[l + l#lni^] 

arc area 

3na^ 

2 ™^ [1 + 1^] 


Table I. Comparison between cycloids on a plane and on a 
sphere. In plane case, the curtate and prolate cycloids are 
traced by a point at radii b < a and h > a, respectively. 

Like a cycloid on plane, in addition to the instanta¬ 
neous eigenstate of the initial Hamiltonian, any initial 
states corresponding to a point on the rim of the rolling 
circle generates the cycloid. Initial states inside and out¬ 
side of the rolling circle trace the curtate cycloid and the 
prolate cycloid on the Bloch sphere, respectively. The 
arc length and arc area of the cycloid on the sphere are 
obtained by Bjelicapii. Table U shows the comparison of 
cycloids on a plane and on a sphere. As shown by the 
blue curve in Fig. [U the curve traced by the axis of the 
rolling circle is parallel to the base line, i.e., the adiabatic 
path. The axis path can be interpreted as an example 
of transitionless drivingi^^— to accelerate the adiabatic 
evolution. This could be understood by considering the 
axis path as a new evolution path of the spin, which 


is driven non-adiabatically by the parallel-rotating mag¬ 
netic field, not by a slowly rotating magnetic field along 
the axis path. Precisely, this can be done by adding the 
transitionless-driving Hamiltonian Ho{t) which cancels 
the non-adiabatic effect of the original time-dependent 
Hamiltonian Hit). So the new time-dependent Hamilto¬ 
nian H {t) -|- Ho {t) drives non-adiabatically a quantum 
state along the adiabatic path of H{t). This technique 
could be used in flipping neutrons non-adiabatically by 
arrange flipping In order to drive the spin-up state (north 
pole) to the spin-down state (south pole) along the merid¬ 
ian line 

Non-adiabatic resonance. The non-adiabatic term, the 
second part in Hamiltonian Q, causes the quantum tra¬ 
jectory deviate from the adiabatic path. The question 
we address now is how the quantum evolution follows 
the adiabatic path as the rotation of the magnetic field 
becomes slower. Even before approaching the adiabatic 
limit, an evolved state could end up the adiabatic target 
state under some condition. This is so called adiabatic 
resonance. Let us consider the magnetic field is rotated 
by angle (3 in time T, i.e., ujT = /3. Clearly, the Bloch 
vector r(t), which is initially aligned to the z-axis in the 
adiabatic frame (that is, the instantaneous eigenstate at 
t = 0), will point again to the z-axis if a Bloch vector in 
the adiabatic frame is completely rotated n times, i.e., 
^ = 27m. The non-adiabatic resonance condition is 
given by 


where n is the number of cycloid arcs. 

Arc area of a cycloid as geometric phase Let us take a close 
look at physics associated with the geometric properties 
of a cycloid curve on a Bloch sphere. The length and area 
are two basic geometric quantities of a curve. First, we 
discuss the area enclosed by the cycloid curve. To this 
end, we consider again the magnetic field which rotates 
completely around z-axis by azimuthal angle 6 as shown 
in Fig. [TJ As is well known, when the magnetic field ro¬ 
tates slowly enough, the evolved state remains in the in¬ 
stantaneous eigenstate, Eq. o, and accumulates, in ad¬ 
dition to the dynamical phase. Berry phase which is pro¬ 
portional to the solid angle enclosed by the adiabatic path 
n(t).— When the adiabatic condition is released, a cycled 
quantum state accumulates the dynamical phase and the 
geometric phase, called Aharonov-Anandan (AA) phase, 
which is proportional to the area of the curve of the quan¬ 
tum evolution in the projected Hilbert spacei^. When 
the resonance condition meets, the Bloch vector v{t) re¬ 
turns to its initial position after a complete rotation of 
the magnetic field, so the cycloid curve is closed, we 
calculate the AA phase and explore how the AA phase 
approaches to Berry phase in adiabatic limit. With the 
exact solution ©, AA phase is written as 

lAA = TT -I- 7d (9a) 
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where 


(AE)T - I 


{AE)T ~ 


Id = 


h 


dt = 


— — COS a 
A 


(9b) 


In the adiabatic limit, i.e., A —>■ 0, the dynamic phase jd 
becomes jd —t —7rcos0. Thus the AA phase converges 
to the Berry phase 


lim 7 Aa[A] = 7Berry = 7r(l - COS 6») . 


( 10 ) 


As shown in Fig. [U the difference between A A and Berry 
phases is the sum of the arc areas of a cycloid. It is 
inversely proportional to running time T so AA phase 
becomes Berry phase in adiabatic limit as expected. 

Length of a cycloid as time-energy uncertainty and its singu¬ 
lar limit Let us turn to physics related with the length of 
the cycloid curve. The Heisenberg position-momentum 
uncertainty relation is based on the anti-commutation re¬ 
lation between two operators. However, the energy-time 
uncertainty is different because time in quantum mechan¬ 
ics is not an operator. Anandan and Aharonov^i gave a 
nice interpretation of the energy-time uncertainty rela¬ 
tion as the distance of the quantum evolution measured 
by the Fubini-Study metric in the projective Hilbert 
space. The length L of the quantum evolution between 
two orthogonal states, here from |0) to |1), is expressed 
as 


L = 2 




(11) 


where AE = is the uncertainty 

in energy during running time T. The length of any curve 
connecting two orthogonal states |0) and |1) is greater 
than or equal to the shortest distance between them, i.e., 
the geodesic line of the length tt, L > tt. So the energy¬ 
time uncertainty is given by 


{AE)T > 


h 


( 12 ) 


This tells that the minimum time T required for trans¬ 
forming to an orthogonal state is bounded by T > 
h/4:{AE), as shown by Mandelstam and TamniS^, Flem- 
ing2^, VaidmanSij and Levitin and Toffoli^^. 

Now let us examine how the length of the quantum evo¬ 
lution changes as the speed of a rotating magnetic field 
from the north to the south poles along the geodesic line 
varies, as shown in Fig. [21 The adiabatic theorem dictates 
that if the magnetic field rotates slowly enough, the quan¬ 
tum evolution is well approximated by an instantaneous 
eigenstate. As depicted in Fig. |2l in adiabatic limit, the 
path of the quantum evolution approaches to the adia¬ 
batic path, i.e., the geodesic line. One would expect that 
the length of the quantum evolution in adiabatic limit 
becomes just that of the adiabatic path (geodesic line) 
because the difference between the two paths, measured 
by the enclosed area (difference between AA and Berry 



Figure 2. Trajectories, infidelity, and length. Top panel shows 
two trajectories in red on a Bloch sphere for 1/2A = 1 and 
7, respectively. The blue line is the adiabatic path or the 
trajectory of a magnetic field. Middle and bottom panels 
plot the infidelity, the probability deviating from |1), and the 
distance of the quantum evolution, respectively, as a function 
of the adiabatic parameter A. Blue points in bottom panel 
represent for the perfect transition. 


phases), becomes vanished. So the time-energy uncer¬ 
tainty in adiabatic limit would be minimum. However, 
that is not the case. In the adiabatic limit the length 
of the quantum evolution becomes L = A, but not tt as 
explained below. This is called the diagonal paradox^^ 
or the limit paradoxal in calculus. Some well-known ex¬ 
amples showing singular limits^^ are the classical limit 
of quantum mechanics, and the limit of zero viscosity^ 
called d’Alembert’s paradox. 

With the exact solution of the Schrddinger equation, 
one can calculate the length of quantum evolution. After 
some algebra, the length L as a function of adiabatic 
parameter A is given by 


m = 


yrTA^ 


77 v^1 + a2 



r 

X \/ 

1 - 

Jo V 



(13) 


cos^ X + sin^ X 


dx 


In the limit of A —>■ 0, while the integrand of Eq. (fT51) be¬ 
comes smaller, the interval of integration become larger. 
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So one obtains 


lim L[A] = 4 (14) 

A^O 

which is greater than the geodesic length tt between the 
north and the south poles. This result can be also under¬ 
stood in terms of the product of the arc length of a cycloid 
and the number of cycloid arcs needed. At non-adiabatic 
resonance condition with D,T = 27rn and uiT = tt, the 
radius a of the cycloid is given by a = ^ = ^. 

In adiabatic limit, the cycloid can be seen as a plane 
cycloid, so the arc length is just 8a. Thus the length 
becomes L ^ 8a x n = A. 

Fig. m plots the infidelity, the probability of deviating 
from the target state |1), and the distance L of quan¬ 
tum evolution as function of the adiabatic parameter 
1/A = ujojuj. At the first non-adiabatic resonance, i.e., 
1/2A = 1, the curve is composed of a single cycloid and 
is deviated from the adiabatic path. However, its length 
is lightly bigger than tt. This implies {AE)T m h/A. On 
the other hand, in adiabatic limit, i.e., 1/2A ^ 1, the 
curve is composed of many cycloids with smaller radius 
and get close to the adiabatic path with length tt. In 
adiabatic limit of A —?► 0, while the red curve in Fig. [2] 
is approaching to the blue one, the adiabatic path, its 
length converges to 4 not tt, as shown before. Then the 
time-energy uncertainty becomes {AE)T ~ h/A which 
is not minimum. The quantum adiabatic theorem tells 
that the instantaneous eigenstate (|4|) is a good approxi¬ 
mation to the true quantum evolution if the Hamiltonian 
changes slowly enough. While the non-adiabatic geomet¬ 
ric phase converges to the adiabatic geometric phase, the 
time-energy uncertainty doesn’t. This shows that the 
adiabatic limit is singular in the sense that the instan¬ 
taneous eigenstate cannot capture all physical properties 
in that limit. 

Discussion 

In conclusion we have shown that the Bloch vector of 
a spin in a rotating magnetic field traces a cycloid on a 
Bloch sphere. Like on plane, different initial states trace 
prolate and curtate cycloids, and the trajectory parallel 
to the adiabatic path. The perfect non-adiabatic reso¬ 
nance is geometrically interpreted as the complete rolling 
of a cycloid. Two fundamental geometric quantities, the 
area enclosed by a cycloid curve and its length, are con¬ 
nected to two physical quantities, the geometric phase 
and the time-energy uncertainty, respectively. The arch 
areas of a cycloid gives rise to the difference between 
AA phase and Berry The energy-time uncertainty be¬ 
comes {AE)T Tn hi TT which is greater than the minimum 
energy-time uncertainty h/A. We found the quantum adi¬ 
abatic limit is singular, similar to the diagonal paradox 
or d’Amlembert’s paradox in the limit of zero viscosity. 
In adiabatic limit, while the A A phase converges to Berry 
phase, the length, time-energy uncertainty, doesn’t con¬ 
verge to that of the adiabatic path. 


In mathematics, the isoperimetric inequality tells the 
relation between the circumference L of a closed curve 
and the area A it encloses. The isoperimetric inequality 
on a sphere is obtained by Levy2£ and is given by 

L^>A{Att-A), (15) 

where the equality holds if and only if the curve is a 
circle. In our case, it tells the relation between the ge¬ 
ometric phase and the time-energy uncertainty. For a 
slow rotation of the magnetic field along the great circle 
passing the north to the south poles, the area is given 
by A = 27r, and the quantum state acquires the phase 
factor, —1. If the quantum evolution were the great cir¬ 
cle, its length would be L = 27r according to the equality 
condition of isoperimetric inequality. The distance of the 
quantum evolution is 8, but not 27r, although it looks like 
a great circle. 

The geometric phases, the area enclosed by a cycloid 
curve, have been measured with various spin-1/2 systems 
such as a neutron in a rotating magnetic field, the polar¬ 
ization of light in a coiled optical fiber, and qubits. The 
time-energy uncertainty, the length of a cycloid curve, 
and its singular limit can be measured with same sys¬ 
tems. While the geometric phase is measured via inter¬ 
ference between an evolved and initial quantum states, 
a trajectory on a Bloch sphere seems to be needed for 
calculating the time-energy uncertainty because of diffi¬ 
culty in measuring energy fluctuation. However, with the 
rapid advancement in manipulating a qubit, it is possible 
to track a trajectory of a qubit on a Bloch sphere. Espe¬ 
cially, we notice that Roushan et al^ traced the cycloid 
curve on a Bloch sphere in the experiment of measuring 
the non-adiabatic geometric phase with superconducting 
qubits. 

Adiabatic approximation is one of fundamental theo¬ 
rem in quantum mechanics, so it has many applications, 
for example, Born-Oppenheimer approximation and adi¬ 
abatic quantum computing. The results here could give 
an opportunity to deepen our understanding of adiabatic 
approximation. 


Methods 

Time-evolution of a spin in a rotating magnetic fields. 

The evolved quantum state at time t is given by 

m)) = AHt)u{t)Amm). 


where U(t) is the time-evolution operator in the adiabatic 
frame. For the magnetic field rotating about the z axis by 
angle 6, the effective Hamiltonian in the adiabatic frame 
is given by 


HeS 


hwQ hw -ay 

cr^ - — [cos 0a ^ - sint'cTj;) 

hfi 

-m • (T 

2 


where H = \/uJq + 2u}u;o cos6 + m = cosai — 
sina£ = TUzZ -|- nixX is the direction of the effective 
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magnetic field in the adiabatic frame, = cosa = {ip{t)\H'^{t)\ip{t)) = and 

(wq + w cos 0)/fl, and rux = — sina = —to sind/Q. Thus 

the time-evolution in the adiabatic frame is given by ((V'(OI^(^)IV'(O)) 

g—i/Tefft/fi _ gi^y-m-CT 

Calculation of time-energy uncertainty The magnetic held 


3/4 


COS 


ftt 


I / 2 2 \ • 2 

-f [m^ — m^) sm 


is rotated from the north pole to south pole along 
the geodesic line. The initial state is |0). With 
the exact solution, it is straightforward to calcu¬ 
late the length L = j-AE{t)dt with AE{t) = 

- {{tp{t)\H{t)\tp{t))f. One obtains 


where ml —ml. = (wg—a;^)/(wQ-|-w^) = (1 —A^)/(l-|-A^). 
By changing the variable in the integrand, one obtains 
Equation (USD. 

Calculation of time-energy uncertainty The trajectory of 
the Bloch vector (^/!(t)|cr|'!/i(t)) is plotted with the exact 
solution and by solving the time-dependent Schrddinger 
equation numerically with the Runge-Kutta method. 
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